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Preface 

This  thesis  is  an  investigation  of  the  present  capabilities  of  a 
specific  Space  Shuttle  Upper  Stage  Vehicle.  It  is  also  an  attempt  to 
make  a contribution  towards  the  economical  use  of  the  vehicle  in 
terms  of  various  types  of  energy  management.  The  approach  to  the 
problem  is  general  in  nature  so  that  it  may  be  extended  to  a wide  range 
of  shuttle  missions  and  vehicles. 

1 would  like  to  acknowledge  those  who  have  helped  and  supported 
me  throughout  this  project.  My  thanks  go  to  my  advisor.  Professor 
Gerald  M.  Anderson,  for  his  aid,  especially  in  the  initial  phase  of  the 
study.  1 am  also  indebted  to  Dr.  Richard  M.  Potter  of  the  Electrical 
Engineering  Department  whose  recommendations  and  interest  kept  me 
going  through  difficult  times. 

Most  of  the  credit  for  this  modest  work  must  be  given  to  my  wife, 
Suzy.  Her  love  and  devotion,  along  with  personal  sacrifices,  are  what 
made  this  final  product  possible.  It  is  to  Suzy,  and  our  three  children, 
Suzanne,  Kara,  and  John  that  this  thesis  is  dedicated. 


i 

; 

i 


Contents 


Page 


Preface ii 

List  of  Figures v 

List  of  Tables vi 

List  of  Notations  vii 

Abstract xi 

I.  Introduction 1 

Background  1 

Statement  of  the  Problem 2 

Assumptions 4 

II.  Formulation  of  Multipoint  Boundary  Value  Problems  ...  7 

Coplanar  Boundary  Value  Problems 8 

Case  One  12 

Case  Tvjo 12 

Case  Three 13 

Noncoplanar  Boundary  Value  Problem 13 

lU.  Numerical  Algorithm 21 


Background  

Algorithm 

ORBTRAN  . . . 

NSOIA 

CALFUN  .... 
DFEQ  ..... 

SLOPE  

Program  Integration 


IV.  Discussion  of  Results 31 

Results  of  Coplanar  Transfers 32 


iii 


«*■ 


Contents 


Page 

Thrust  Termination 

Effects  of  Varying  Transfer  Angle 

Effect  of  Offloading  Fuel 

Effects  of  Changing  Terminal  Orbit  Parameters  . . 

General  Comments 

Results  of  the  Noncoplanar  Transfers 

Cross  Check  of  Noncoplanar  and  Coplanar 

Transfers 

Effect  of  Changing  Inclination  Angle 

Use  of  Integration  Scheme  to  Check  Noncoplanar 
Problem 

V.  Conclusions 55 

Bibliography 58 

Appendix  A;  Coplanar  Equations  of  Motion  and  Boundary 

Conditions 59 

Appendix  B:  Noncoplanar  Equations  of  Motion  and  Boundary 


Conditions 55 

Appendix  C:  Derivation  of  Necessary  Conditions 73 

Appendix  D:  Computer  Listing  of  Coplanar  Problem 80 

Appendix  E:  Computer  Listing  of  Noncoplanar  Problem  ....  88 

Vita 97 


List  of  Figures 


Page 


1-1  The  Space  Shuttle  Vehicle 


1 


3- la  Numerical  Algorithm  Flow  Chart 29 

3-  lb  Numerical  Algorithm  Flow  Chart 30 

4- 1  Effect  on  Payload  of  Varying  Transfer  Angle 

(ef  = 0) 38 

4-2  Effect  on  Payload  of  Varying  Transfer  Angle 

(ef  0) 39 

4-3  Effect  on  Payload  of  Varying  Transfer  Angle 

(Pf  = 14,500  nm)  41 

4-4  Effect  on  Payload  of  Offloading  (Geosynchronous 

Mission) 43 

4-5  Effect  on  Payload  of  Offloading  (12  Hour  Mission)  . 45 

4-6  Time  History  of  Control  (Geosynchronous 

Mission) 48 

4-7  Time  History  of  Control  (12  Hour  Mission) 49 

A-1  Coplanar  Thrust  Vector 60 

B-1  Noncoplanar  Transfer 65 

B-2  Noncoplanar  Thrust  Vector 66 

B-3  Spherical  Coordinate  System 67 

B-4  Spherical  Angles  for  Noncoplanar  Problem 70 


List  of  Tables 


Table  Page 

I Boeing  Burner  II  Specifications 5 

II  Summary  of  Coplanar  MPBVPs 14 

III  Summary  of  Noncoplanar  MPBVPs . 20 

IVa  Results  of  Thrust  Termination  Problems  33 

IVb  Results  of  Thrust  Termination  Problems 34 

V Varying  p£  of  Terminal  Circular  Orbit  (0=180°) 45 

VI  Maximum  Payload  Transfers  to  Elliptical  Orbits 45 

VII  Coplanar  Solutions  for  Orbital  Transfers  (Case  Three).  . 50 

VIII  Payload  Check  of  Coplanar  and  Noncoplanar  (i=0) 

Problems  53 

DC  Effect  of  Inclination  Angle  on  Payload  ( At£  = 120.  89).  ...  53 

X Results  of  Integration  Scheme  (i=10,  0*=14O,  ^tj  = 139.  52)  . 55 


vi 


List  of  Notation 


English  Symbols 


convergence  criterion 


AFIT 


Air  Force  Institute  of  Technology 


effective  exhaust  velocity 


CMU 


canonical  mass  unit 


CTU 


canonical  thrust  unit 


differential  change 
degrees 


step  length  used  in  search  routine 


canonical  distance  unit 


eccentricity 


unit  vector 


error  function  vector 


force  due  to  thrust 


force  due  to  gravity 

specific  angular  momentum  vector 


hamiltonian 


inclination  angle  between  two  orbits 


specific  impulse 


vii 


lUS 


Interim  Upper  Stage 


T 


.1 


J 

J 

k 


h 


Ibf 

Ibw 

m 

MPBVP 

nm 

P 

r 


SAMSO 


sec 

SRM 

SSV 

STS 

t 

h 

^2 


i 

cost  function  I 

Jacobian  matrix 
switching  function 

direction  cosine  of  thrust  vector  i=l,  2,  3 
pounds -force 

pounds -weight  | 

vehicle  mass  j 

multipoint  boundary  value  problem  j 

nautical  mile 

semi-latus  rectum 

radial  distance  from  center  of  earth 

radius  at  apogee 

radius  at  perigee 

radius  of  earth 

Space  and  Missile  Systems  Organization 
seconds 

solid-rocket  motor 
Space  Shuttle  Vehicle 
Space  Transportation  System 
time 

time  of  first  stage  cutoff 
time  of  second  stage  ignition 
time  of  staging 


viii 


T 


TPBVP 

TU 


u 

USV 


X 

X 


thrust  magnitude 

two-point  boundary  value  problem 
canonical  time  unit 
control  variable 
upper  stage  vehicle 

velocity  components  in  "e^, "e^,  directions 
respectively 

state  variable  vector 

vector  of  unknowns 


Greek  Symbols 


0 

7 

5 


Ml 

^^Imax 

At£ 

^*-2max 

e 

e' 

\ 


1? 


thrust  steering  angle 

spherical  angle 

spherical  angle 

a small  variation 

elapsed  time  for  first  stage 'burn 

maximum  allowable  time  for  first  stage  burn 

elapsed  time  for  second  stage  burn 

maximum  allowable  time  for  first  stage  burn 

central  transfer  angle  in  initial  orbital  plane 
measured  counterclockwise  from  launch 

central  transfer  angle  in  final  orbital  plane 
measured  counterclockwise  from  launch 

costate  vector 

gravitational  parameter 

lagrange  multiplier  vector  associated  with  terminal 
conditions 


IX 


lagrange  multiplier  vector  associated  with 
constraints 

spherical  angle  measured  clockwize  from  Z 


Subscripts 

actual 

d 

£ 

o 


computed  value 
desired  value 
final  value 
initial  value 
vector  quantity 


I 

! 


Superscripts  j 

— unit  vector 

• first  time  derivative 

• • second  time  derivative 

~ augmented  function 

• optimal  quantity 


I 


X 


GA/MC/76D-11 


i= 


Abstract 

I 

A two- stage,  solid-fuel  proposal  for  the  Interim  Upper  Stage 
Vehicle  of  the  Space  Transportation  System  has  been  selected  for  use 
by  the  Department  of  Defense.  It  is  the  purpose  of  this  study  to 
investigate  the  capabilities  of  such  a vehicle  in  terms  of  the  maximum 
payload  which  can  be  delivered  to  orbit.  This  optimal  payload  prob- 
lem is  examined  in  light  of  three  energy  management  techniques.  The 
first  technique,  thrust  termination,  involves  shutting  off  engine  thrust 
prior  to  complete  use  of  propellant.  The  second  technique  investi- 
gates the  effects  on  payload  of  varying  the  central  angle  through  which 
the  transfer  is  made.  The  third  technique,  offloading,  examines  the 
possibility  of  reducing  the  amou.it  of  available  fuel  for  either  stage 
prior  to  the  mission  to  determine  if  payloads  can  be  increased.  Finite 
burn  periods  are  assumed  in  this  study.  A multipoint  boundary  value 
problem  results  from  the  appearance  of  interior  constraints  in  the 


problem.  A numerical  technique  is  used  to  generate  solutions  for  a 
range  of  transfers  in  the  special  coplanar  problem.  Conclusions  are 
discussed  for  each  of  the  three  energy  management  techniques  in  the 
coplanar  transfer.  Even  though  non-optimal  noncoplanar  results  were 
obtained  it  is  believed  that  the  general  coplanar  results  can  be  extended 
to  the  noncoplanar  problem. 
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MAXIMUM  PAYLOAD  ORBITAL  TRANSFERS 
FOR  A SPACE  SHUTTLE  SOLID-FUEL 
UPPER  STAGE  VEHICLE 


I.  Introduction 


Background 

The  Space  Transportation  System  (STS)  is  a concept  which  is 
designed  to  reduce  the  costs  of  earth  orbital  operations  and  advance 
the  capabilities  of  the  United  States  in  space.  The  STS  consists 
of  a Space  Shuttle  Vehicle  (SSV),  an  upper  stage  vehicle  (USV)  and  the 
necessary  support  systems.  The  SSV,  depicted  in  Fig.  1-1,  is  designed 
to  deliver  scientific  and  defense  satellites  to  earth  orbits. 


Fig.  1-1.  The  Space  Shuttle  Vehicle 
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The  orbiter  vehicle,  shown  in  Fig.  1-1,  is  a reusable  component 
of  the  SSV  which  is  delivered  to  a low-Earth  orbit  through  the  use  of 
an  expendable  external  liquid-fuel  tank  and  two  reusable  solid  rocket 
boosters.  The  SSV  has  the  capability  to  deliver  the  orbiter  to  low 
earth  orbits  of  altitudes  up  to  500  nautical  miles  (nm).  However,  many 
of  the  missions  for  which  the  STS  is  designed  require  higher  orbits. 

For  these  missions  a propulsive  upper  stage  vehicle  (LiSV),  which  is 
carried  in  the  cargo  bay  of  the  orbiter,  is  deployed  by  the  orbiter  in 
the  low  earth  orbit.  The  USV  is  then  used  to  deliver  payloads  beyond 
the  capability  of  the  orbiter  vehicle.  The  STS  was  designed  to  relax 
the  weight  and  volume  constraints  which  previously  have  been  placed 
on  payloads  transported  to  high  earth  orbits. 

The  USV  that  is  presently  being  developed  is  the  Interim  Upper 
Stage  (lUS).  The  lUS  vehicle  is  capable  of  delivering  single  or  multiple 
payloads  to  high  earth  orbits  by  use  of  a single  lUS  vehicle.  Boeing’s 
Burner  II  version  of  the  lUS,  a two-stage,  solid-fuel  vehicle,  has  been 
selected  by  the  Space  and  Missile  Systems  Organization  (SAMSO)  for 
use  in  the  initial  phases  of  the  Shuttle  program. 

Statement  of  the  Problem 

One  way  to  examine  the  economical  use  of  the  lUS  for  an  orbital 
transfer  mission  is  to  determine  the  maximLun  payload  that  can  be 
delivered  to  orbit.  This  thesis  addresses  this  maximum  payload  prob- 
lem for  one  specific  lUS  vehicle.  Three  energy  management  techniques 
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are  used  to  examine  this  optimization  problem  in  an  effort  to  study  I 

the  capabilities  of  the  lUS.  The  objective  of  this  optimization  problem 
is  to  determine  the  time  history  of  the  controls  to  give  the  maximum 
payload  transfers. 

The  first  technique,  early  thrust  termination,  implies  incom- 
plete use  of  the  available  propellant.  This  can  be  achieved  by  blowing 

i 

off  the  engine  nozzle  or  reducing  the  pressure  in  the  combustion  ^ 

chamber  of  an  engine.  The  two-stage  lUS  proposals  submitted  to 

SAMSO  do  not  have  this  capability,  but  the  concept  is  examined  here  | 

i 

to  determine  if  discarding  excess  fuel  in  orbit  would  enable  more  pay- 

load  to  be  delivered.  I 

1 

The  second  energy  management  scheme  is  offloading,  that  is,  ' 

reducing  the  amount  of  propellant  for  a given  stage  prior  to  the  mis- 
sion. Examining  the  optimization  problem  with  the  added  dimension 
of  offloading  provides  a sizing  analysis  in  terms  of  the  maximum  pay- 
load  that  can  be  delivered  to  orbit  for  different  amounts  of  fuel.  Thus 
for  a specific  payload,  the  "best"  selection  of  available  propellant 
loadings  is  possible. 

The  third  energy  management  technique,  varying  the  transfer 
angle,  concedes  the  fact  that  there  is  excess  fuel  which  must  be  carried 
and  used.  The  question  now  becomes  one  of  how  to  use  that  excess 
fuel.  The  present  concept  is  to  use  a reaction  control  system  to  burn 
extra  propellant.  However,  the  concept  of  varying  the  transfer  angle 
is  presented  here  as  an  alternative  to  such  an  additional  control  system 
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which  is  designed  to  burn  excess  fuel. 

The  orbital  transfers  made  in  this  study  use  finite  periods  of 
thrusting  rather  than  the  more  common  impulsive  approximation  of 
velocity  changes.  A finite  burn  is  viewed  as  a more  exact  solution  to 
the  problem  as  it  considers  the  fact  that  thrust  direction  may  not  be 
constant  during  thrusting  periods.  Another  factor  resulting  in  a more 
realistic  solution  for  the  finite  burns  is  that  the  finite  burn  case  con- 
siders the  gravitational  forces  which  must  be  overcome  during  periods 
of  thrusting.  The  maximum  payloads  for  the  finite  burn  problem  would 
thus  be  expected  to  be  somewhat  lower  than  those  for  the  impulsive 
approximations. 


Assumptions 

Upon  consideration  of  both  solid-  and  liquid-fuel  lUS  proposals, 
SAMSO  has  granted  a contract  to  Boeing  Aerospace  Company  for  its 
Burner  II  solid-rocket  motor  (SRM)  version  of  the  lUS.  This  partic- 
ular two- stage  vehicle  does  not  have  the  capability  to  thrust  terminate, 
so  that  concept,  as  developed  in  this  study,  pertains  to  future  applica- 
tions of  such  a vehicle.  The  vehicle,  however,  does  presently  have 
the  capability  of  offloading  fuel.  The  specifications  of  the  Burner  II 
vehicle  can  be  found  in  Table  I.  It  is  assumed  that  the  SRMs  provide 
a constant  thrust  magnitude  with  a maximum  amount  of  burning  time 
for  each  stage. 
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Table  I 

Boeing  Burner  II  Specifications 


r 


Stage 

1 Summary 

Thrust 

41,  757.4 

Ibf 

Specific  Impulse 

291.  8 

sec 

Structure  Weight 

2,437.  0 

Ibw 

Propellant  Weights 

17,  300.  0 

Ibw 

or  19,967.0 

Ibw 

Stage 

2 Summary 

Thrust 

14, 076. 6 

Ibf 

Specific  Impulse 

300.  8 

sec 

Structure  Weight 

1,  362. 0 

Ibw 

Propellant  Weight 

4,  700. 0 

Ibw 

The  controls  which  we  use  to  optimize  the  problem  are  the  thrust 
directions.  For  the  thrust  termination  problems,  the  numerical 
algorithm  discussed  in  Chapter  II  also  seeks  the  optimal  times  for 
first  stage  cutoff  and  second  stage  ignition.  The  equations  which  gov- 
ern the  motion  of  the  vehicle  are  those  as  derived  in  an  inverse-square 
gravitational  field.  Accelerations  from  sources  other  than  the  SRMs 
and  the  gravitational  force  of  the  earth  are  assumed  to  be  negligible. 
The  assumption  is  also  made  that  changes  in  body  attitude  do  not  affect 
the  motion  of  the  vehicle. 

For  the  orbital  problems  studied  that  involve  transfers  to  final 
elliptical  orbits,  the  major  axis  of  the  final  orbit  is  assumed  to  be 
aligned  with  the  point  on  the  initial  orbit  at  which  the  transfer  begins. 
This  assumption  allows  the  true  anomaly  in  the  polar  equation  of  a 


} 
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conic  orbit,  Eq  (A-16),  to  be  equated  to  the  central  angle  through  which 
the  transfer  is  made,  6 . 

The  low  earth  parking  orbit  of  the  orbiter  vehicle  is  assumed 
throughout  this  thesis  to  be  a circular  orbit  around  the  earth  at  an 
altitude  of  150  nm. 


II.  Formulation  of  Multipoint  Boundar 


Value  Problems 


The  objective  of  this  study  is  to  determine  the  time  histories  of 
the  controls  for  a nonlinear  system  in  such  a manner  that  the  initial 
mass  is  maximized  subject  to  boundary  constraints.  The  problem  of 
thrust  termination  mentioned  in  Chapter  I is  atypical  in  that  boundary 
constraints  in  this  problem  appear  not  only  at  initial  and  terminal 
points,  but  also  at  interior,  or  intermediate  points.  The  initial  con- 
straints are  given  by  the  parameters  of  a parking  orbit  from  which  the 
mission  begins.  The  terminal  constraints  are  specified  a priori  and 
they  define  the  orbit  to  which  the  transfer  is  made.  The  intermediate 
constraints  pertain  to  physical  parameters  of  the  vehicle.  Two  of 
these  interior  constraints  are  placed  on  the  maximum  amount  of  time 
that  each  motor  of  the  two-stage  lUS  can  burn.  A third  interior  con- 
straint involves  the  mass  which  the  vehicle  jettisons  at  the  staging 
point.  These  three  inequality  constraints  are  discussed  in  Appendix  C. 
They  necessitate  the  formulation  of  a multipoint  boundary  value  prob- 
lem (MPBVP)  because  of  the  possibility  of  discontinuities  in  the  vari- 
ables or  equations  at  interior  points. 

It  is  also  shown  in  Appendix  C that  three  separate  cases  can  be 
considered  for  the  thrust  termination  problem.  The  first  two  cases 


are  concerned  with  terminating  thrust  of  one  of  the  two  lUS  stages 
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before  the  propellant  loaded  for  that  stage  has  been  conapletely  used. 
The  third  case  implies  complete  use  of  propellant  in  both  stages.  The 
problem  of  offloading  fuel  is  covered  by  case  three  for  which  a differ- 
ent amount  of  propellant  is  assumed  to  be  available  for  that  stage. 

The  multi-point  boundary  value  problem  for  each  of  the  three 
cases  that  result  from  applying  the  necessary  conditions  of  Appendix 
C are  presented  here. 


Coplanar  Boundary  Value  Problems 

The  state  equations  which  govern  the  coplanar  orbital  transfer 
are  derived  in  Appendix  A and  are  as  follows: 


(2-1) 


li  = Xr  cos  O' 
m 


C 


where 


r is  the  distance  from  the  center  of  the  earth 
6 is  the  angle  from  the  initial  point  to  the  present  position 
Vj^  is  radial  velocity 


h is  specific  angular  momentum 
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H is  the  gravitational  parameter 


i 

t 

i 

i 
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T is  the  magnitude  of  thrust 

m is  the  vehicle's  mass  at  a given  time 

a is  the  thrust  steering  angle  measured  clockwise 
from  the  local  horizontal 

C is  the  effective  exhaust  velocity 


Since  we  want  to  maximize  the  payload  for  a fixed  amount  of 
propellant,  the  cost  function,  J,  is: 


J = -m(tQ) 


(2-2) 


where  tQ  represents  initial  time.  The  mathematics  of  Appendix  C 
seeks  a minimum  to  this  cost  function.  The  minimum  of  a negative 
initial  mass,  with  a fixed  change  in  mass,  is  therefore  the  same  as 
maximizing  final  mass.  The  Hamiltonian,  H,  is  defined  as: 


H = Xj,Vr  + X,  Jl  + X, 


R 


r h^  _ + X sin  «1 

L m J 


(2-3) 


+ Xu  ~ cos  a - X X 
“ m C 


where  ^V|^*  ^h*  ^m  costate  variables. 

Application  of  the  optimality  condition  of  Appendix  C,  Eq  (C-21), 
and  enforcement  of  the  strengthened  Legendre  condition,  yields  the 
following  expression  for  the  thrust  steering  angles: 
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r 


sin  Of  = 


(2-4) 


where  P = w ^ 

1 ^ R h 


The  costate  equations  for  this  problem  are  defined  by  Eq  (C-20). 
Substituting  the  controls  of  Eq  (2-4)  yields: 


±_ 

r3  . 


2 A^h  + 


X,  2rT 


A,  = 0 


^h  - 3 


i - -TP 

“ 2 
m 


(2-5) 


An  initial  circular  parking  orbit  at  an  altitude  above  the  earth 
of  150  nm  is  assumed  throughout  this  study.  The  initial  conditions 
are  therefore: 


r«  = R + 150  nm 

V 0 

<>0  = 0 


(2-6) 
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(2-6) 


Vr  = 0 


where  R is  the  radius  of  the  earth  and  ( ) denotes  value  at  the 
w o 

initial  point.  Another  initial  condition  is  given  for  the  mass  costate 
by  Eq  (C-  17)  as: 


'm.  = 1 


(2-7) 


The  terminal  constraints  are  repeated  here  from  the  derivation 


in  Appendix  C: 


1 + e£  cos  0£ 


^£  - specified 


e£h£  sin  0f 


(2-8) 


h£  = fjl^ 


where  ( )£  denotes  the  value  at  the  final  time,  p is  the  semi-latus 
rectum  and  e is  eccentricity.  We  also  know  from  the  transversality 
conditions  derived  in  Appendix  C that  the  mass  costate  at  the  final 
time  is  given  by: 


(2-9) 
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The  remainder  of  the  boundary  value  problem  must  be  devel- 


I 


I 


oped  for  each  of  the  three  cases  described  in  Appendix  C. 

Case  One.  This  case  is  characterized  by  incomplete  use  of  the 
second  stage  fuel  while  the  first  stage  is  completely  burned.  The 
necessary  conditions  unique  to  this  problem  are  described  by 
Eq  (C-11).  From  these  conditions  we  know  the  value  of  the  Hamil- 
tonian at  the  final  time  is: 


I 

'!■ 


r 

!■ 

r 

i 


H(t£)  = 0 (2-10) 

It  is  also  known  that  the  Hamiltonian  is  continuous  at  second  stage 
ignition.  The  partial  derivative  of  the  Hamiltonian  with  respect  to 
thrust  at  that  time  must  therefore  be  zero.  The  expression  described 
by  this  relationship  is  known  as  the  switching  function.  The  equa- 
tion is: 


i 


(2-11) 


where  t^  is  the  time  of  second  stage  ignition. 

This  MPBVP  is  summarized  in  Table  II.  It  can  be  seen  that  a 
solution  to  this  problem  with  twelve  equations  and  twelve  unknowns 
is  possible. 

Case  Two.  This  case  involves  thrust  termination  of  only  the 
first  stage  of  the  lUS.  The  necessary  conditions  derived  in  Appendix 
C,  Eq  (C-13)  yields  the  following  condition: 


H(tf)  = H4(t2)  - H3(t2) 


(2-12) 


where  H3(t2)  is  the  Hamiltonian  just  prior  to  second  stage  ignition  and 
H4(t2)  is  the  Hamiltonian  an  instant  after  ignition  of  the  second  stage. 

The  Hamiltonian  being  continuous  at  the  first  stage  cutoff,  tj,  yields: 

The  boundary  value  problem  is  summarized  in  Table  U. 

t 

Case  Three.  The  third  case,  which  serves  as  a basis  of  compar- 
ison for  the  other  two  cases,  involves  the  complete  burning  of  all 
loaded  propellant  in  the  lUS.  The  comparison  with  the  case  three 
problem  is  used  to  indicate  if  thrust  termination  of  either  stage  yields 
any  extra  payload.  This  case  is  also  used  for  the  analysis  of  offload- 
ing fuel.  The  associated  necessary  condition  derived  in  Appendix  C is: 

H(tf)  = H4(t2)-H3(t2)  (2-14) 

Since  the  Hamiltoniein  is  discontinuous  at  both  tj  and  t2,  there  is  no  ( 

i 

equation  involving  the  switching  function.  The  MPBVP  is  given  in  j 

Table  II.  As  in  the  other  two  cases,  it  can  be  seen  that  the  number 
of  unknowns  match  the  number  of  equations. 

Noncoplanar  Boundary  Value  Problem  | 

The  general  problem  of  a noncoplanar  orbital  transfer  is  | 

depicted  in  Fig.  B-1.  The  appropriate  state  equations  for  this  prob-  ■ 

! 

lem  as  derived  in  Appendix  B are:  ] 
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f = V 


R 


V 


R 


sin^<^  - 


m sin<f> 


V.  - -VrV^,  +V^  Sind,  T^3 

V<^  - -!- COS  <f)  + 


m 


(2-15) 


where  r is  the  radius  from  the  earth's  center  to  the  lUS 

0 is  the  transfer  angle  measured  in  the  initial  orbital  plane 
<f,  is  the  angle  measured  from  the  vertical  axis  to  the  vehicle. 
Vr,  Vg,  and  are  the  velocity  state  variables 
m is  the  vehicle's  mass  at  a given  time 
M is  the  gravitational  parameter 
T is  thrust  magnitude 

and  are  the  direction  cosines  of  the  thrust  vector 
With  the  cost  function,  J,  for  this  problem  defined  as 

J=-m(tp)  (2-16) 


the  Hamiltonian  becomes: 
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where  V associated  with 

the  state  variables. 

As  in  the  coplanar  problem,  the  optimality  necessary  condition 
and  the  strengthened  Legendre-Clebsch  condition  applied  to  the 
Hamiltonian  yield  the  following  expressions  for  the  thrust  directions; 

P 

= “Ayg- . (2-18) 

P sin^ 

f,  - 

P 

where  P= 

^ y sin4>/ 

The  necessary  condition  for  the  costates,  defined  by  Eq  (C-20),  is 
applied  to  the  above  Hamiltonian  with  the  expressions  for  thrust 
direction  substituted  in  the  proper  places.  The  results  are; 
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" + " ^V^('^R^</>“V  sin<^cos(/))] 

'^e  = 0 

2 

A|^=  sin<^cos^-2Ay^  V^V^  - T cos  ^ {cos^4>  - 

r 5 : — 7 — ; — 

r sin^<#>  m p sin-^  <#>  r 

sin^</>) 

\ =i  " ^R  (2-20) 

R j;*  ^ 

- — [x«+  2 Ay  V,9  sin^  <j>  - Ay  Vj^  - ZXy  cot  </>  + 2'-y,  V^  sin<;>cos<^] 

7p</>-  c°t</.  + Av^Vj^] 

TP 


The  assumption  of  an  initial  circular  parking  orbit  once  again 
defines  some  of  the  initial  conditions  as  shown  below.  The  terminal 
conditions  derived  in  Appendix  C are  also  shown  here. 


^o  ” ^ 


<t>o  = 90° 

= 0 ’ (2-21) 
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= _21L  sin  i sin  0f 

I »• 


> (2-22) 


where  ( indicates  value  at  initial  point 
( )£  denotes  value  at  terminal  point 

fl£  is  the  transfer  angle  in  the  plane  of  the  terminal  orbit 
h is  angular  momentum 
e is  eccentricity 

i is  the  inclination  angle  between  the  initial  and  final  orbits 
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Additional  boundary  conditions  on  the  mass  costate  are  derived 


in  Appendix  C.  Those  boundary  conditions  are: 

(to)  = 1 (2-23) 

A„,(tf)  = 0 (2-24) 

As  in  the  coplanar  problem,  a specific  set  of  necessary  equa- 
tions exists  for  each  of  the  three  classes  of  problems  studied  here. 
The  reader  is  referred  to  the  formulation  of  the  coplanar  MPBVP 
earlier  in  the  chapter  for  an  explanation  of  those  additional  necessary 
conditions.  A summarization  of  the  MPBVP  for  the  nonplanar  case 
is  found  in  Table  III. 
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Table  III 

ary  of  Noncoplanar 


Total 


III.  Numerical  Algorithm 


Background 

Having  derived  the  necessary  conditions  for  an  optimal  trajec- 
tory and  formulated  a MPBVP,  a method  to  find  a solution  for  the 
problem  must  now  be  generated.  The  necessary  conditions  given  in 
Chapter  II  must  be  satisfied  at  each  point  in  the  time  interval  of 
interest.  A closed-form  solution  of  this  problem  could  not  be  found. 
Therefore,  a numerical  technique  is  used  to  search  for  the  unknowns 
such  that  the  necessary  equations  and  boundary  conditions  are  satis- 
fied. That  which  follows  is  an  explanation  of  the  numerical  algorithm 
developed  to  solve  the  MPBVP. 

Algorithm 

Due  to  the  possibility  of  discontinuities  in  the  problem  vari- 
ables, the  MPBVP  formulated  in  Chapter  II  is  separated  into  three 
related  parts  that  can  be  implemented  in  a computer.  The  first  part 
covers  the  first  thrusting  period  from  the  initial  point  to  first  stage 
cutoff  at  time  tj.  The  next  portion  covers  the  coasting  period  of  the 
trajectory  from  first  stage  cutoff  until  second  stage  ignition.  It  is 
during  this  time  interval  that  staging  occurs.  The  last  of  the  three 
parts  includes  the  second  stage  burning  period.  The  endpoint  of  this 
problem  signifies  the  final  time,  t£,  at  which  point  the  terminal 
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boundary  conditions  must  be  satisfied. 

The  numerical  algorithm  formulated  in  this  thesis  consists  of 
several  routines  which  perform  unique  functions  within  the  main 
algorithm.  An  explanation  of  the  processing  involved  in  those  rou- 
tines and  the  program  integration  are  included  in  this  chapter. 

ORBTRAN.  This  routine  in  the  main  program  sets  up  a nor- 
malization routine,  accepts  input  data,  and  forms  an  initial  guess  for 
the  unknowns.  ORBTRAN  also  initiates  the  calls  to  the  subroutines 
in  the  program. 

The  normalization  routine  used  in  this  study  is  similar  to  the 

scheme  proposed  in  Ref  3.  All  distance,  velocity,  time,  mass  and 

thrust  units  in  this  problem  are  divided  by  canonical  units  prior  to 

their  use  in  the  search  routine.  The  distance  unit,  DU,  represents 

the  radius  of  the  earth.  A circular  reference  orbit  is  assumed  just 

above  the  surface  of  the  earth.  Velocity  and  time  canonical  units  are 

therefore  defined  by  the  speed  of  a hypothetical  satellite  in  that  refer- 

3 

ence  orbit.  The  gravitational  parameter,  n,  is  therefore  1 

Tu2 

Consistent  mass  and  thrust  units  are  also  developed.  The  mass  unit, 
CMU,  that  provided  the  best  convergence  was  5000  Ib-sec^/ft.  The 
thrust  unit,  CTU,  was  defined  as  follows 

IDU 

1 CTU  = 1 CMU  X T (3-1) 

ITU'^ 

The  next  main  function  of  ORBTRAN  is  to  receive  inputs  in  the 
form  of  vehicle  parameters,  initial  conditions,  and  terminal 
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conditions.  Specific  impulses  of  each  stage  are  input  to  ORBTRAN. 
These  values  are  used  to  determine  the  effective  exhaust  velocities 


i 

I 

i 

I 


for  the  motors. 

The  initial  conditions  for  the  problem  are  determined  from  the 
parking  orbit  parameters  which  are  input  to  ORBTRAN.  Those 
parameters  are  eccentricity  and  semi-latus  rectum.  Terminal  con> 
ditions  are  determined  in  ORBTRAN  by  inputs  of  the  transfer  angle, 
semi-latus  rectum  and  eccentricity  of  the  desired  terminal  orbit. 

The  equations  for  the  terminal  conditions  listed  in  the  Appendices  are 
used  to  determine  the  values  the  states  should  have  at  t£.  All  bound- 
ary conditions  are  then  normalized  within  ORBTRAN. 

Initial  guesses  of  the  unknowns  are  input  to  the  main  program. 
Once  the  correct  values  of  the  unknowns  have  been  determined  for  a 
particular  mission  those  values  are  then  used  as  inputs  for  another 
similar  orbital  transfer.  This  recursive  use  of  the  algorithm  pro- 
vides good  guesses  for  the  unknowns  and  thus  significantly  reduces 
the  required  computational  time.  ORBTRAN  receives  as  outputs 
from  NSOIA  the  values  of  the  unknowns  that  yield  a solution  to  the 
MPBVP. 

NSOIA.  This  subroutine,  which  is  on  permanent  file  at  the 
AFIT  computer,  was  developed  by  M.  J.  D.  Powell.  (See  Ref  7). 
NSOIA  receives  as  inputs  a vector  of  unknowns,  x,  and  several 
parameters  which  direct  its  operation.  The  algorithm  was  developed 


to  numerically  solve  a system  of  n non-linear  equations  in  n 
unknowns : 

ik  (xi,X2,X3  . . . = 0 k=l,...,n  (3-2) 

where_f  is  a vector  of  error  functions  (non-linear  equations)  in  the 
unknown  parameters,  x^.  It  is  claimed  in  Ref  7 that  this  algorithm 
maintains  a fast  convergence  rate  in  addition  to  the  capability  to 
achieve  convergence  if  initial  guesses  of  the  unknowns  are  poor 
estimates. 

NSOIA  calculates  the  Jacobian  matrix,  J^,  of  the  error  function: 


A Newton-Raphson  step  gives: 

f(x)+J5x  = 0 (3-4) 

where  8x  is  a correction  to  x.  It  is  therefore  possible  to  find  a Sx 
for  a given  step  length  so  that  the  following  equation  is  met: 

f(x  + 6x)  C _f(x)  (3-5) 

A correction  vector  is  found  at  each  iteration  until  the  desired  accur- 
acy is  achieved,  that  is: 

n 

E [1k(x)]^<acc 
K=1 
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For  this  thesis,  with  the  variables  normalized  to  canonical  units, 

- 12 

ACC  was  specified  as  10 

Another  important  parameter  in  NSOIA  which  must  be  properly 
chosen  to  insure  convergence  is  the  increment  size  for  the  unknowns, 
DSTEP,  A step  length  that  is  too  small  would  require  an  excessive 
number  of  iterations  to  achieve  convergence.  On  the  other  hand, 
DSTEP  must  be  small  enough  to  insure  that  Eq  (3-5)  is  satisfied  at 
each  iteration.  A value  of  10  for  DSTEP  proved  to  be  the  best 
choice  for  the  problems  of  this  thesis. 

The  unknowns  in  the  problem,  x,  are  all  incremented  at  the 
same  constant  DSTEP  within  NSOIA.  The  proper  scaling  of  those 
unknowns  is  therefore  required.  Improper  scaling  could  result  in  a 
negligible  change  in  the  error  functions  at  each  increment  in  the 
unknowns,  (see  Eq  (3-3)).  Such  a result  would  make  convergence  to 
the  proper  solution  nearly  impossible.  Eor  example,  it  was  deter- 
mined early  in  the  study  that  the  MPBVP  was  very  sensitive  to 

changes  in  the  initial  costates.  They  had  to  be  incremented  with  a 
n 

DSTEP  of  10”  . Such  a change  in  the  unknown  initial  mass  parameter 
at  each  iteration  would  yield  negligible  changes  in  the  error  functions 
due  to  that  parameter  since  it  has  a value  in  the  range  of  30,  000  Ibw. 
Convergence  to  a solution  would  be  virtually  impossible  in  this  case. 
The  normalization  scheme  mentioned  earlier  places  the  proper 
significance  on  the  increment  size  for  the  unknowns  at  each  iteration 
and  thus  eliminates  some  of  the  convergence  difficulties. 


1 
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NSOIA  internally  calls  subroutine  CALFUN  which  is  used  to 
find  the  error  functions.  A more  detailed  description  of  NSOIA  can 


be  found  in  Ref  7, 

CALFUN.  The  main  task  of  this  subroutine  is  to  provide  an 
error  function,  f,  for  use  in  NSOIA.  The  components  of  the  error 
function  are  determined  from  the  difference  between  the  desired  and 
computed  values  of  the  equations  for  a specific  case  of  the  MPBVP. 
For  example  for  case  three  of  the  coplanar  problem  we  have 

^ ^d  ^actual 

where  Br  and  Br  are  the  desired  and  computed  values  for  the 

d ^actual 

transfer  angle. 

As  inputs  CALFUN  receives  the  desired  values  of  the  boundary 
conditions,  the  initial  guesses  on  the  unknowns  and  the  vehicle 
parameters.  CALFUN  then  acts  in  an  executive  manner  to  integrate 
the  differential  state  and  costate  equations  for  each  portion  of  the 
trajectory.  The  actual  values  of  the  states  which  are  required  for 
the  error  functions  are  then  extracted  from  the  integration  routine 
at  the  proper  points. 

The  integration  routine  in  CALFUN  includes  three  different 
calls  to  the  integration  scheme.  The  first  call  is  related  to  the  first 
thrusting  period.  First  stage  thrust  magnitude,  effective  exhaust 
velocity  and  tj  are  required  for  this  routine.  The  switching  function, 
k (t2)  is  determined  at  the  end  of  the  first  portion  of  the  MPBVP. 
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The  next  portion  of  the  MPBVP,  representing  the  coast  period, 
receives  the  final  values  of  the  states  and  costates  from  the  first 
part  as  starting  values.  Thrust  and  effective  exhaust  velocity  are 
both  zero  during  coasting.  Upon  initialization  of  these  differential 
equations,  the  first  stage  structure  weight  is  subtracted  from  the 
mass  state  variable.  At  the  end  of  this  period  at  t7,  the  switching 
function,  k (t^),  and  Hamiltonian,  H3(t2)  are  determined. 

The  final  values  from  the  coasting  period  are  used  as  inputs 
for  the  states  and  costates  in  the  third  portion  of  the  MPBVP. 

The  Hamiltonian,  H4(t2),  is  determined  after  this  portion  of 
the  MPBVP  is  initialized.  Since  this  phase  covers  second  stage 
thrusting,  the  appropriate  thrust  and  effective  exhaust  velocity  must 
be  provided.  The  required  terminal  boundary  conditions,  ^m{t£), 
and  H(tf)  are  determined  at  the  endpoint  of  this  problem. 

DFEQ.  This  routine  is  also  on  permanent  file  at  the  AFIT 
computer.  It  is  comprised  of  two  subroutines,  SET  and  STEP. 

These  subroutines  accomplish  the  integration  of  the  first-order 
differential  equations  in  the  problem.  The  inputs  required  by  these 
subroutines  are  the  initial  values,  total  integration  time,  and  step 
size.  SET  is  first  called  to  initialize  the  integration  scheme.  Then 
STEP  is  called  to  integrate  the  equations  of  motion  through  a time 
interval,  DT.  A Runge-Kutta  method  is  used  to  integrate  over  three 
steps  after  the  initial  point  is  determined  by  SET.  A fourth-order 
Adams-Bashforth-Adams-Moulton  predictor-corrector  scheme  is 
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then  applied  to  the  next  two  points.  This  procedure  uses  the  equa- 
tions from  subroutine  SLOPE  to  integrate  over  the  entire  interval,  t. 

SLOPE.  Subroutine  slope  sets  up  the  differential  state  and 
costate  equations  as  derived  in  Chapter  11.  The  equations  are  gen- 
eral in  nature  so  they  may  be  used  for  either  thrust  or  coast  periods 
of  a trajectory.  The  inputs  to  this  subroutine  are  the  values  of  the 
states  and  costates  at  a given  time  when  called  from  DFEQ.  This 
subroutine  computes  the  slope  of  a given  state  or  costate  which  is 
subsequently  used  by  DFEQ  to  estimate  the  value  of  that  parameter 
at  some  future  time,  t + 4t. 

Program  Integration 

The  overall  integration  of  these  subroutines  to  comprise  the 
numerical  algorithm  can  be  seen  in  Fig.  3-1.  The  main  program, 
ORBTRAN,  performs  the  executive  function  of  the  algorithm.  It 
provides  the  call  to  the  search  routine,  NSOIA.  That  routine  then 
changes  the  initial  guesses  of  the  unknowns  such  that  the  error  func- 
tion decreases  in  magnitude.  Internal  to  NSOIA  is  a call  to  CALFUN 
which  directs  the  integration  routine  so  that  the  error  function  can  be 
determined.  NSOIA  repeatedly  calls  CALFUN  until  the  error  func- 
tion is  within  the  desired  accuracy.  The  outputs  of  NSOIA  are  the 
unknowns  which  yield  the  desired  accuracy  in  the  error  function. 


ORBTRAN 


Fig.  3- la.  Numerical  Algorithm  Flow  Chart 


Fig.  3- lb.  Numerical  Algorithm  Flow  Chart 


IV*  Discussion  of  Results 
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Convergence  to  a solution  to  the  MPBVP  proved  to  be  very 
sensitive  to  changes  in  the  initial  estimates  of  the  unknowns.  Those 
unknowns  include  the  initial  values  of  the  costates  and  the  unknown 
times.  The  algorithm  failed  to  converge  to  a solution  for  poor 
guesses  of  those  unknowns.  The  differences  between  the  desired  and 
computed  values  of  the  terminal  conditions  were  not  within  acceptable 
limits  when  poor  initial  estimates  of  the  unknowns  were  used.  Since 
physical  reasoning  could  not  be  used  to  obtain  guesses  for  the  A 
unknowns,  a great  deal  of  computer  time  was  required  to  achieve 
convergence  for  the  initial  problem.  Once  a solution  was  obtained, 
however,  the  correct  values  of  the  unknowns  were  used  as  inputs  for 
a similar  mission  with  a small  change  in  the  transfer  angle  or  in  the 
radius.  In  this  manner  a range  of  solutions  were  obtained. 

Both  the  coplanar  and  noncoplanar  results  are  presented  in  this 
chapter.  Analytical  justifications  and  validation  of  results  are  also 
included.  The  maximum  payload  orbital  transfer  results  are  pre- 
sented for  the  two  primary  design  reference  missions  of  the  lUS, 

The  first  of  these  is  a transfer  to  a circular  geosynchronous  orbit 
at  an  altitude  of  19,  323  nm.  The  second  mission,  known  as  the  12 
hour  mission,  defines  the  terminal  orbit  as  an  elliptical  orbit  with 
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an  apogee  altitude  of  21,450  nm  and  a perigee  altitude  of  350  nm. 
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The  results  are  also  presented  for  the  two  main  first  stage  propellant 
loadings  listed  in  Table  I.  The  first  loading  of  19,  967  pounds  gives 
a maximum  burn  time  of  139.  52  seconds,  while  the  17,  300  pound 
loading  enables  120.89  seconds  of  first  stage  burn. 

Results  of  Coplanar  Transfers 

Results  for  each  of  the  three  energy  management  techniques 
were  obtained  for  the  coplanar  transfer  problem. 

Thrust  Termination.  At  the  beginning  of  this  study,  it  was 
believed  that  the  ability  to  thrust  terminate  might  be  desirable  in 
terms  of  increasing  payload  capability.  However,  successful  orbital 
transfers  were  not  obtained  for  either  the  case  one  or  case  two 
problem.  Numerous  approaches  were  attempted  to  obtain  a solution 
to  the  thrust  termination  cases.  The  missions  were  changed  to 
determine  if  thrust  termination  could  produce  solutions  for  only  a 
certain  class  of  missions.  Different  sets  of  the  initial  estimates  of 
the  unknown  costates  were  used  as  inputs,  but  to  no  avail.  The  best 
accuracies  that  could  be  achieved  for  either  case  of  thrust  termina- 
tion were  obtained  by  using  the  values  of  the  unknowns  for  a success- 
ful case  three  transfer  in  which  all  the  fuel  was  used.  Those  results 
appear  in  Table  IV.  It  is  apparent  from  an  examination  of  the  burning 
times  and  initial  masses  listed  in  the  table  for  the  case  one  and  case 
two  problems  that,  even  when  case  three  results  are  input  to  the 
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Table  IVa 

Results  of  Thrust  Termination  Problems 
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thrust  termination  problems,  the  algorithm  does  not  always  drive 
towards  the  case  three  solution.  For  instance,  the  case  one  problem 
for  a 12  hour  mission  with  B — 180®  shows  a second  stage  burn  of 
about  52  seconds  which  is  far  below  that  burning  constraint.  The 
established  convergence  criterion  was  never  achieved  for  the  thrust 
termination  problems.  If  early  thrust  ternnination  is  not  optimal, 
it  appears  that  a solution  should  at  least  have  been  obtained  by  the 
algorithm  driving  to  the  limiting  case  three  solution.  The  reason 
for  the  failure  of  the  algorithm  to  provide  such  solutions  most  prob- 
ably lies  in  the  fact  that  case  one  and  case  two  have  error  functions 
which  include  switching  functions.  The  sensitivity  of  these  functions 
to  the  costate  values  at  a particular  time  is  very  great,  as  discovered 
in  Ref  6.  The  error  function  in  the  algorithm  involving  the  switching 
function  was  the  primary  contribution  to  the  relatively  high  value  of 
final  accuracy  which  was  on  the  order  of  10”®  for  the  case  one  prob- 
lem. For  the  problems  involving  those  additional  necessary  condi- 
tions it  was  therefore  not  possible  to  achieve  the  desired  accuracy. 

A variety  of  missions  and  transfer  angles  were  tested  with  the 
results  being  similar  to  those  in  Table  IV.  Both  of  the  first  stage 
propellant  loadings  were  also  tested.  Convergence  to  a solution 
failed  in  every  instance. 

Although  convergence  to  a solution  was  not  obtained  for  the 
thrust  termination  problems,  an  analytical  argument  is  presented 
here  to  support  the  fact  that  thrust  termination  cannot  increase 
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payload  capability.  Solutions  to  the  case  three  problem  are  avail- 
able for  fuel  loadings  which  are  less  than  the  maximum  loadings. 
Those  results  can  be  used  to  determine  if  case  one  or  case  two  can 
increase  payload.  For  example,  with  18,  900  lbs  of  first  stage  fuel 
and  4,  700  lbs  of  second  stage  fuel  (maximum  second  stage  loading), 
the  case  three  payload  capability  is  about  5,  100  lbs.  If  the  second 
stage  loading  is  decreased  to  4,400  lbs  for  the  same  first  stage 
loading,  the  case  three  problem  indicates  that  the  final  mass  is 
increased  to  5,  300  lbs.  In  the  case  one  analysis  of  this  latter  prob- 
lem, however,  the  burning  of  4,  400  lbs  of  fuel  is  interpreted  as 
incomplete  use  of  available  second  stage  propellant.  This  means 
that  300  lbs  of  the  final  mass  is  unused  fuel,  or  deadweight.  The 
payload  is  therefore  5,000  lbs.  This  indicates  a decrease  in  payload 
when  all  available  fuel  is  not  used.  This  result,  which  consistently 
holds  in  such  an  analysis  for  decreasing  second  stage  burn,  shows 
that  terminating  second  stage  thrusting  early  will  not  increase  pay- 
load  capability.  The  optimal  solution  to  the  case  one  problem,  which 
has  the  possibility  to  use  less  than  all  the  loaded  second  stage  fuel, 
is  always  to  use  the  maximum  amount  of  available  fuel. 

Similar  results  are  obtained  in  the  analysis  of  the  case  two 
problem.  Terminating  the  first  burn  early  decreases  the  payload 
capability.  For  example,  with  the  maximum  first  and  second  stage 
fuel  being  used  approximately  5,  900  lbs  of  payload  can  be  delivered. 
Keeping  the  second  stage  burn  time  constant  and  terminating  first 


stage  thrusting  7.  52  seconds  early  decreases  the  final  mass  to  5,  100 
lbs.  Subtracting  the  unused  first  stage  fuel  from  the  final  mass  yields 
a payload  of  4,  000  lbs.  The  conclusion  that  can  be  arrived  at  is  that 
the  optimal  solution  to  the  case  two  problem  is  to  burn  the  maximum 
amount  of  time  on  the  first  burn. 

Effects  of  Varying  Transfer  Angle.  In  contrast  to  the  thrust 
termination  problems,  a variety  of  solutions  were  obtained  for  the 
third  case  of  the  multi-point  boundary  value  problem.  Thus  the 
other  two  energy  management  techniques  of  varying  transfer  angle 
and  off-loading  fuel  can  be  examined.  One  cross-check  of  the  accur- 
acy of  this  problem  was  obtained  by  a comparison  with  the  Hohmann 
transfer  problem.  For  both  first  stage  propellant  loadings  the 
Hohmann  approximation  for  the  initial  mass  was  within  one  percent 
of  the  maximum  initial  mass  for  the  finite  burn  solution  from  the 
computer. 

The  results  of  this  portion  of  the  study  do  not  present  any 
startling  new  ideas.  The  findings  of  this  aspect  of  the  problem  agree 
with  the  classical  Hohmann  transfer  problem  which  was  developed  in 
the  1920's.  However,  a few  interesting  observations  can  be  made  for 
the  finite  burn  problem  and  its  possible  applications  for  an  lUS  mis- 
sion. First  of  all,  the  transfer  angle  for  the  finite  burn  problem  at 
which  we  have  maximum  payload  is  greater  than  180°.  This  fact  is 
evidenced  in  Figs.  4-1  and  4-2.  This  result  is  due  to  the  fact  that  the 
transfer  angle  changes  a small  amount  during  the  finite  burn.  The 
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Fig.  4-2.  Effect  on  Payload  of  Varying  Transfer  Angle  (e£  0) 
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angle  change  through  the  coast  period  remains  approximately  180^^. 
Maximum  payload  transfers  therefore  occur  at  about  185°.  This 
result  is  consistent  throughout  the  range  of  transfers  examined. 

The  results  of  varying  the  transfer  angle  seem  to  imply  that 
this  technique  may  be  used  as  an  energy  management  scheme.  For 
a mission  with  a payload  less  than  that  which  can  be  delivered  w ith 
a given  propellant  loading,  the  present  lUS  concept  is  to  use  a reac- 
tion control  system  to  burn  excess  fuel.  An  examination  of  Figs. 

4-1  and  4-2  show  that  such  an  additional  control  system  may  be  unnec- 
essary in  its  application  to  burn  excess  fuel.  Changing  the  transfer 
angle  of  the  mission  can  achieve  the  same  objective  of  burning  all  of 
the  available  propellant  in  an  orbital  transfer  mission. 

It  is  also  interesting  to  examine  the  range  of  transfer  angles 
for  a given  mission.  A comparison  of  the  graphs  in  Fig.  4-1  shows 
that  an  increase  in  first  stage  propellant  loading  not  only  increases 
the  payload  which  can  be  delivered,  but  it  also  increases  the  range 
of  transfer  angles.  The  range  of  transfer  angles  to  the  elliptical 
orbit  shown  in  Fig.  4-2  is  smaller  than  that  for  transfer  to  geosynch- 
ronous orbit.  Figure  4-3  shows  the  range  of  transfer  angles  for  a 
mission  to  a 14,500  nm  orbit.  In  this  case  not  only  has  payload  capa- 
bility increased,  but  the  range  of  transfers  has  also  increased  signif- 
icantly. 

Effect  of  Offloading  Fuel.  It  was  mentioned  in  Chapter  I that 
Boeing's  Burner  II  version  of  the  lUS  has  the  capability  to  offload 
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fuel.  The  failure  of  the  thrust  termination  problem  to  produce  any 
positive  results  brought  forth  the  question  of  offloading.  It  was 
originally  thought  that  if  thrust  termination  failed  to  increase  payload 
that  the  case  three  problem  with  less  propellant  available  should  also 
enable  less  payload  to  be  delivered  to  orbit.  The  results  of  the 
investigation  appear  here  in  the  form  of  a sizing  iinalysis.  For  a 
mission  with  a specified  payload  it  is  possible  to  examine  such  an 
analysis  to  determine  how  much  propellant  is  required  to  achieve  the 
terminal  orbit.  Economic  factors,  engineering  considerations  and 
other  factors  can  then  be  considered  to  determine  which  sizing  is 
best  for  a given  mission.  The  sizing  analysis  was  performed  for  the 
two  primary  missions  of  the  lUS  with  a transfer  angle  of  180^.  The 
resvilts  can  be  found  in  Figs.  4-4  and  4-5. 

The  results  of  the  sizing  analysis  for  the  geosynchronous  mis- 
sion show  that  offloading  first  stage  propellant  does  not  enable  more 
payload  to  be  transported.  A decrease  in  first  stage  burn  will  also 
decrease  the  available  payload.  The  results  for  offloading  second 
stage  propellant  prove  to  be  a bit  more  interesting.  When  the  maxi- 
mum available  propellant  is  loaded  for  the  first  stage,  offloading 
second  stage  fuel  does  not  yield  any  extra  payload  to  orbit.  However, 
as  shown  in  Fig.  4-4,  as  the  first  stage  fuel  is  decreased  offloading 
second  stage  propellant  is  advantageous  in  terms  of  payload  capability. 
With  a first  stage  burn  of  132  seconds  the  maximum  payload  is 
delivered  with  a second  burn  of  97.5  seconds.  Decreasing  the  second 
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Fig.  4-4.  Effect  on  Payload  of  Offloading  (Geosynchronous  Mission) 


burn  another  1.  5 seconds  (total  offloaded  fuel  of  206  pounds)  would 
still  enable  more  payload  to  be  delivered  than  that  possible  with  the 
full  second  stage  loading.  At  the  lower  limit  of  the  first  stage  pro- 
pellant loadings  (iitj  = 120.  89  seconds)  the  maximum  payload  capa- 
bility occurs  at  4t^  = 87.6  seconds.  Offloading  as  much  as  1029 
pounds  of  second  stage  fuel  for  that  first  stage  loading  w’ould  still 
enable  the  same  payload  to  be  delivered  as  that  at  full  second  stage 
loading. 

The  sizing  analysis  for  the  12  hour  mission  reveals  significantly 
different  findings.  The  offloading  of  first  stage  propellant  still  does 
not  yield  any  extra  payload.  However,  offloading  second  stage  fuel 
increases  payload  for  einy  first  stage  loading.  It  might  be  mentioned 
here  that  the  transfer  trajectories  for  the  orbital  transfer  through  an 
angle  of  180°  are  elliptical  orbits  with  relatively  high  values  for 
eccentricity.  It  would  appear  that  the  velocity  change  required  at  the 
second  burn  would  not  be  very  large  for  insertion  into  a similar 
elliptical  orbit.  This  is  exactly  the  result  that  was  obtained.  As 
shown  in  Fig.  4-5,  the  maximum  payload  is  delivered  in  the  neighbor- 
hood of  a 30  second  burn  for  the  second  stage  motor.  For  these  12 
hour  missions  70  percent  of  the  available  fuel  can  be  offloaded  and 
still  deliver  more  than  3,  000  pounds  of  additional  payload  to  orbit. 
Offloading  another  700  pounds  of  fuel  would  still  enable  the  same 
amount  of  payload  to  be  transported  to  this  elliptical  orbit. 
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Effects  of  Changing  Terminal  Orbit  Parameters.  In  addition 


f. 

ij 

V to  the  analysis  of  the  three  previously  mentioned  energy  management 

^‘i 

techniques,  the  maximum  payloads  that  can  be  delivered  to  different 
terminal  orbits  is  also  investigated.  The  semi-latus  rectum,  p£,  of 

I the  final  circular  orbit  was  varied  and  the  results  appear  in  Table  V. 

1 

! The  expected  result  is  that  more  payload  can  be  delivered  to  an  orbit 

' ! with  a smaller  radius. 


Table  V 

Varying  pf  of  Terminal  Circular  Orbit  (®=  180°) 


pf  (nm) 

5000 

10000 

13000 

16000 

19000 

19500 

21000 

Payload  (lbs) 

8611 

6329 

6048 

5356 

4471 

4350 

4023 

Transfers  were  also  made  to  elliptical  terminal  orbits.  Chang- 
ing the  apogee  distance,  r^,  and  perigee  distsince,  r^,  of  the  elliptical 
orbit  is  described  in  this  problem  by  eccentricity,  ef,  and  semi-latus 
rectum,  pf.  The  results  for  the  transfers  to  elliptical  orbits  for  trans- 
fer angles  of  180°  are  presented  in  Table  VI. 

Table  VI 

Maximum  Payload  Transfers  to  Elliptical  Orbits 


j 

(: 

j' 

I 


Terminal 

Orbit 

Parameters 

«f 

Pf  (nm) 
r^  (nm) 
IV.  (nm) 

.218 

19,489. 

20,000. 

19.000. 

.0517 

20,685. 

22,000. 

19.500. 

. 1525 
18,965. 

23.000. 

16. 000. 

.2445 

15,778. 

23.000. 

12.000. 

.456 
10, 325. 
22,000. 
6.000. 

.593 
8,088. 
24,  893. 
3.793. 

Payload  (lbs) 

5,597 

5,192 

5,017 

5,  193 

5,218 

4,724 

General  Comments.  The  objective  of  this  optimization  as 
described  in  Chapter  I is  to  determine  the  time  histories  of  the  con- 
trols for  the  transfer  trajectory  which  yields  maximum  payload.  The 
mathematical  relationships  for  the  controls  are  given  by  the  optimal- 
ity condition.  Graphs  of  the  time  histories  of  the  controls  for  the 
two  primary  lUS  missions  are  given  in  Figs.  4-6  and  4-7. 

Due  to  the  sensitivity  of  the  numerical  algorithm  to  the  unknown 
lagrange  multipliers  a great  deal  of  computer  time  was  required  to 
obtain  solutions  for  the  coplanar  problems.  The  recursive  use  of  the 
algorithm  helped  reduce  the  amount  of  computer  time,  but  converg- 
ence to  solutions  was  still  costly.  In  the  transfer  to  the  geosynchro- 
nous orbit  for  example,  the  full  range  of  transfers  required  well  over 
2000  seconds  of  execution  time.  This  time  does  not  include  the  time 
required  for  convergence  to  the  initial  transfers  angle  in  that  range. 
Table  VII  shows  how  the  unknowns  in  the  problem  varied  for  the 
different  missions. 

A pure  integration  routine,  using  the  values  of  the  unknowns 
from  a successful  orbital  transfer,  was  used  as  a check  on  the  accur- 
acy of  the  numerical  algorithm  to  search  for  the  correct  unknowns. 
The  results  of  the  integration  routine  show  that  the  desired  final 
conditions  are  met  for  the  orbital  transfer  problems.  Distance 
errors  between  desired  and  computed  values  of  terminal  conditions 
are  on  the  order  of  ten  feet.  Errors  in  angles  are  approximately  10~^ 
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degrees.  Velocity  errors  are  on  the  order  of  10'^  feet  per  second. 
Accuracy  even  better  than  that  obtained  is  possible,  to  a certain 
extent,  by  demanding  more  accuracy  from  the  search  routine, 

NSOIA. 

Results  of  the  Noncoplanar  Transfers 

It  was  noted  in  the  previous  section  that  the  numerical  algo- 
rithm was  very  sensitive  to  errors  in  the  unknowns  and  that  converg- 
ence was  often  difficult  to  obtain.  This  problem  became  even  more 
pronounced  in  the  noncoplanar  transfer  due  to  the  additional  states 
and  lag  range  multipliers. 

Cross  Check  of  Noncoplanar  and  Coplanar  Transfers.  The 
first  task  that  was  attempted  for  the  noncoplanar  transfer  was  a 
check  on  the  claim  of  Appendix  B,  That  claim  was  that  the  nonco- 
planar transfer  state  equations  could  be  used  for  a coplanar  problem 
by  simply  setting  the  noncoplanar  inclination  angle  equal  to  zero. 

The  resulting  comparisons  of  the  use  of  two  different  sets  of  equation 
for  the  same  problem  appear  in  Table  VIII,  It  can  be  seen  that  both 
approaches  give  the  same  maximum  payload  to  within  one  tenth  of  a 
pound.  This  holds  true  for  both  first  stage  propellant  loadings. 

Effect  of  Changing  Inclination  Angle.  Due  to  the  extreme  sens- 
itivity to  variations  in  the  unknowns  for  the  noncoplanar  problem,  the 
recursive  use  of  the  algorithm  was  once  again  required.  The  inclina 
tion  angle  was  incremented  by  one  fourth  of  a degree  for  each 


successive  problem.  Approximately  1000  seconds  of  computer  time 
was  still  required  to  acliieve  the  necessary  convergence.  Table  IX 
shows  how  the  unknowns  changed  by  incrementing  the  inclination 
angle. 

The  trends  shown  in  Table  IX  continue  as  the  inclination  angle 
is  increased.  It  therefore  appears  that  it  would  not  be  possible  to 
deliver  a payload  beyond  a five  degree  inclination  angle.  These 
results  are  not  consistent  with  the  claims  of  Ref  4,  that  payloads  of 
more  than  3000  pounds  can  be  delivered  through  an  inclination  angle 
of  28.  5°.  The  numerical  algorithm  was  also  checked  for  a first 
burn  of  139.  52  seconds.  An  inclination  angle  of  10°  and  transfer 
angle  of  140°  were  used  for  the  transfer.  According  to  the  solutions 
from  the  algorithm  only  414  pounds  of  payload  can  be  delivered.  This 
is  once  again  far  below  the  stated  payload  capability  given  in  Ref  4. 
The  next  step  taken  was  to  change  the  transfer  angle  to  determine  at 
what  angle  the  maximum  payload  could  be  delivered.  Changing  the 
transfer  angle  to  145°  for  i = 10°  reduced  the  payload  to  215  pounds. 
Similar  results  were  obtained  when  the  transfer  angle  was  incre- 
mented for  i = 1°.  This  result  is  in  direct  contradiction  of  the  results 
of  the  Hohmann  problem.  Even  though  a noncoplanar  transfer  is 
being  accomplished  increasing  the  transfer  angle  should  result  in 
increased  payload  capability.  Various  other  initial  guesses  for  the 


unknowns  failed  to  yield  convergence  to  the  noncoplanar  transfer. 
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Table  VtU 

Payload  Check  of  Coplanar  and  Noncoplanar  (i=0)  Problems 


Since  the  results  indicate  a non-optimal  solution  for  the  non- 
coplanar  transfer  an  impulsive  approximation  of  a 1°  plane  change 
■was  examined  to  determine  how  m*uch  the  payload  would  be  affected 
by  such  a change.  The  1°  plane  change  for  a Hohmann  transfer  to 
a geosynchronous  orbit  for  the  vehicle  parameters  of  this  problem 
cost  282  lbs  of  payload.  This  is  far  below  the  decrease  in  payload 
for  the  same  plane  change  as  indicated  by  the  computer  results. 

It  thus  appears  that  NSOIA  has  converged  to  a non-optimal  solution. 

Use  of  Integration  Scheme  to  Check  Noncoplanar  Problem.  In 


order  to  determine  just  where  the  search  routine  was  placing  the 
lUS  for  a payload  that  is  obviously  non- optimum,  an  integration 
scheme  was  employed.  The  method  is  the  same  as  that  used  in  the 
coplanar  problem.  The  unknowns  obtained  from  a successful  solu- 
tion in  terms  of  convergence  accuracy  were  input  to  an  integration 
routine.  The  desired  and  actual  values  of  the  terminal  conditions 
were  then  compared.  The  results,  which  appear  in  Table  X,  indi- 
cate that  the  unknowns  are  correct  in  terms  of  getting  the  vehicle 
to  the  proper  terminal  conditions. 


Table  X 

Results  of  Integration  Scheme  (i  = 10,  O'  = l40,  Atj  = 139.  52) 


State 

Computed 
V alue 

Desired 
V alue 

Error 

(Magnitude) 

r 

6.  610734611704  DU 

6.610734442127  DU 

4.  18  ft. 

e 

139.  5675381° 

139. 5675389° 

8. 0x10"^  deg. 

h 

82. 3557286° 

82. 36572992° 

1.  32xl0-^deg. 

Vr 

5.  84655x10-’7  DU/TU 

0.  00 

. 012968  ft/sec 

. 3899?.47944172 
DU/TU 

. 3899240926573 
DU/TU 

. 018155  ft/sec 

04380487048311 

DU/TU 

-.04380155563068 

DU/TU 

. 08559  ft/sec 

Similar  results  were 


obtained  when  the  integration  scheme  was  used  to 


test  the  accuracy  for  the  problem  with  i = 1 


V.  Conclusions 


Three  energy  management  techniques  have  been  applied  to  the 
maximum  payload  orbital  transfers  problems. 

The  coplanar  results  for  these  energy  management  tecliniques 
indicate  that  optimal  solutions  for  thrust  termination  problems  can- 
not be  found  using  the  numerical  algorithm  of  this  study.  However, 
the  analytic  argument  presented  indicates  that  thrust  termination 
does  not  increase  payload  capability  of  the  lUS. 

On  the  other  hand,  varying  the  transfer  angle  and  offloading 
fuel  are  two  possibilities  for  managing  lUS  propellant.  Variation  of 
the  transfer  angle  can  achieve  the  same  objective  as  a reaction  con- 
trol system  that  is  used  to  burn  excess  fuel.  The  sizing  analysis 
provided  by  investigating  the  offloading  problem  enables  one  to 
choose  the  best  of  a variety  of  propellant  loadings  for  a mission 
with  a specified  payload. 

Although  solutions  were  not  obtained  for  the  noncoplanar  prob- 
lem it  is  assumed  that  the  coplanar  results  regarding  the  energy 
management  techniques  can  be  extended  to  the  noncoplajiar  problem. 
That  is,  of  the  three  energy  management  techniques  investigated, 
only  the  thrust  termination  problem  would  not  yield  positive  results. 


Convergence  to  solutions  was  difficult  with  the  non-linear 
equation  solver,  NSOIA,  due  to  the  sensitivity  of  the  problem  to 
such  factors  as  the  costates  and  switching  functions.  A great  deal 
of  care  in  scaling  the  problem  and  choosing  the  best  step  length  for 
NSOIA  still  left  the  fact  that  a great  deal  of  computer  time  was 
required  to  obtain  solutions. 
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Appendix  ^ 


Coplanar  Equations  of  Motion  and  Boundary  Conditions 

The  coplanar  orbital  transfer  as  formulated  here  is  just  a special 
case  of  a noncoplanar  transfer  for  -which  the  inclination  angle  between 
the  initial  and  final  orbits  is  zero.  The  equations  of  motion  derived  in 
this  appendix  provide  state  equations  which  differ  from  those  of  the  non- 
coplanar problem  and  thus  serve  as  a check  for  one  case  of  the  more 
difficult  noncoplanar  transfer. 

Equations  of  Motion 

The  motion  of  a vehicle  that  is  subject  only  to  gravitational  and 
propulsive  forces  can  be  described  by  the  following  differential  equa- 
tion; 


r = G (r,t) 

- — M(t) 


where 

£ is  the  position  vector  of  the  vehicle 

t represents  time 

G(r,  t)  is  the  force  due  to  gravity 

F(t)  is  the  force  due  to  thrust 

M(t)  is  the  mass  of  the  vehicle 

and  ( ) indicate  differentiation  with  respect  to  time 


(A-1) 
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The  assumptions  of  Chapter  I dictate  that  the  gravitational  attraction 
is  an  inverse-square  central  force,  therefore 


G(r,t)  = 

r3 


(A-2) 


•where  p-  is  the  gravitational  parameter  and  r is  the  magnitude  of  the 
position  vector.  Restricting  thrust  to  lie  within  a plane  enables  one  to 
write  F(t)  in  terms  of  thrust  magnitude,  T,  and  steering  angle,  a,  as 
sho-wn  in  Fig.  A-1. 


The  expression  for  acceleration  of  the  vehicle  in  terms  of  radial, 
and  tangential,  e^  , components  becomes: 

£ = ( - + — cos  a)  "ep  + — sina'e«  (A- 3) 

r'^  m m 


Another  expression  for  acceleration  can  be  obtained  by  twice 
differentiating  the  position  vector,  r^r,  of  the  vehicle  with  respect  to 
time.  The  angular  velocity  of  the  vehicle  needed  for  the  differentiation 
is: 


b) 


(A-4) 


■j 
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where  6 is  the  angle  between  r^(t)  and  £ (t  + At),  and  e^  is  the  unit 
vector  normal  to  the  plane  of  motion.  The  well-known  result  of  this 
operation  is 


.2  _ . _ 

£ = (r-r0  )e^  + (2rO+r0)eg 


(A-5) 


Equating  components  of  the  two  acceleration  equations  yields  the  follow- 
ing equations; 


T - lO  = - ~ 


h T 

— + — sin  a 


Zr  6 + T 9 = 


2 m 


T 

— cos  a 
m 


(A-6) 


Three  states  that  may  be  used  in  this  problem  are  position,  r, 
central  angle,  6,  and  radial  velocity,  Vj^.  Two  of  the  differential 
equations  are  therefore; 


r = V 


R 


(A-7  ) 
(A-8) 


where  is  tangential  velocity.  Substituting  these  equations  plus  the 
time  derivatives  of  Eq  (A-8)  into  Eqs  (A-6)  gives  the  following  relation- 
ships: 

2 

r2  m 


Vr  = 


^ +3!.  s in  a 


(A-9) 
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(A-10) 


V.  = — cos  a - 
^ m r 


The  knowledge  that  angular  momentum  is  constant  during  coast  phases 
of  a trajectory  provides  the  motivation  for  its  use  as  a state.  The 
equation  for  angular  momentum  is: 


h = r X V 


(A-11) 


where  h is  angular  momentum  and  _V  is  the  velocity  vector.  Angular 
momentum  for  coplanar  motion  is  therefore 


h = rVfl 


(A-12) 


Substitution  of  Eq  (A-12)  into  Eq  (A-9)  and  Eq  (A-10)  gives  the  final 
form  for  these  state  equations. 

Since  mass  varies  with  time  during  thrusting  periods,  it  can  be 
described  by  a differential  equation  in  terms  of  thrust  and  effective 
exhaust  velocity,  C. 


T 

m 


(A-13) 


The  state  equations  for  coplanar  transfer  therefore  are: 


r = Vr 
a . h 


V - h'"  . T . 

Vr  - ^ sin  a 

r-’ 


(A-14) 
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h - — 
m 


cos  a 


rii  = 


T 

C 


Boundary  Conditions 

The  formulation  of  the  boundary  value  problem  requires  that  the 
values  of  the  states  at  the  endpoints  be  determine  a priori  or  by  the 
numerical  algorithm.  The  initial  values  of  four  of  the  states  are  spec- 
ified by  the  parameters  of  the  parking  orbit  from  which  the  mission 
begins.  Those  conditions  are: 


r(0)  = r^ 
6{0)  = 0 

h(0)  = V 


(A-15) 


The  initial  value  of  the  mass  state  variable  is  unknown  in  that  the 
algorithm  seeks  the  maximum  initial  mass  for  which  a transfer  can  be 
made.  An  initial  guess  at  this  value  was  determined  from  the  specifi- 
cations of  the  Boeing  Burner  II  vehicle. 

It  is  assumed  that  the  final  orbit  parameters  and  the  orbit  trans- 
fer angle  for  a mission  are  determined  a priori.  The  desired  final 
orbit  is  prescribed  by  the  semi-latus  rectum,  p£,  and  eccentricity,  e^. 
For  a given  transfer  angle,  0^,  the  final  radius  is  given  by: 

Pf 

*’f  “ 1 + 6£  cos  (A-16) 
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The  final  radial  velocity  in  terms  of  the  specified  parameters  is; 


'rf  = ef  sin  6^ 

* Pf 


(A-17) 


Final  angular  momentum  is  given  by; 


hf  =>//‘Pf 


(A-18) 


The  final  mass  is  not  required  for  solution  of  the  MPBVP  as  shown  in 
Chapter  II. 


1 


Appendix  i3 

Noncoplanar  Equations  of  Motion  and  Boundary  Conditions 

The  equations  derived  in  this  appendix  can  be  used  for  either  a 
noncoplanar  or  coplanar  transfer  by  adjusting  the  inclination  angle 
between  the  initial  and  final  orbits.  Much  of  the  work  in  deriving  the 
state  equations  and  boundary  conditions  for  this  problem  can  be  found 
in  References  2 and  6 . The  general  case  of  the  noncoplanar  trans- 
fer can  be  seen  in  Fig.  B-1. 


Equations  of  Motion 

As  in  the  coplanar  case  the  motion  of  the  vehicle  is  determined 
by  Eq  (A-1)  which  is  repeated  here: 


The  gravitational  force,  ^(r,  t)  is  written  as  in  Appendix  A,  but  the 
propulsive  force,  ^(t),  is  written  in  terms  of  the  r-</>-0  spherical 
coordinate  system.  Figure  B-2  shows  the  applied  thrust  vector  F(t) 


in  the  r-<t>-6  frame. 


The  equation  of  motion  in  spherical  coordinates  is  therefore 
written  as: 

r = (-Ji+5li)er  + '^^e<^  + Zlle^  (B-2) 

m m m 

where 

M is  the  gravitational  parameter 
r is  the  magnitude  of  the  radius  vector 
T is  thrust  magnitude 
m is  vehicle  mass 

auid  --^1»^2»'^3  direction  cosines  associated  with  a,  P and  y 

' respectively. 

t 

i 
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(B-5) 


..  -2  2 ■ 2 — 
r^=(f-r0  sin 

+ (2  r </>-  r 0 sin  <pcos  </'  + 

+ (2rSsin<#>  + rfl  sin</>  + 2r0‘#>  cos  <j>) 

Having  found  this  second  equation  for  acceleration  it  is  now  possible 
to  determine  the  state  equations.  Differential  state  equations  in  r,  9, 
and  4>  are  defined  as: 


r = Vr 

5 = ^ (B-6) 


r 


Equating  like  components  of  Eq  (B-2)  and  Eq  (B-5)  yields  the  following 
state  equations  for  Vr,  and  : 

Vr  = sin^i>  + ^ - ii  + Hi 

r r r2  rn 

(B-7) 

V^  = sin<l>cos4>  + 

^ r m 

V _ 2V,,V^cos<(>  T13 

r sin<#>  m sin 

An  additional  differential  equation  for  the  time  rate  of  change  of  mass 
is: 

T 


(B-8) 


•where  C is  the  effective  exhaust  velocity. 


The  seven  state  equations  required  for  the  noncoplanar  transfer 
are  given  by  Eqs  (B-6),  (B-7),  and  (B-8). 

Boundary  Conditions 

The  values  of  the  states  at  the  boundaries  must  be  given  or  the 
algorithm  must  seek  their  value  to  give  a solution  to  the  boundary  value 
problem.  Initial  values  •which  are  prescribed  as  in  Appendix  A are 

r(0)  = r^ 

0(0)  = 0 (B-9) 

Vr(0)  = 0 

The  initial  orbit  is  assumed  to  be  in  the  “ex,  ey  plane  as  sho'wn  in  Fig. 
B-2.  The  angle  <f>  is  therefore  given  by 

<;.{0)=90°  (B-10) 

This  assumption  gives  the  value  of  V in  the  initial  orbit  as: 

V^(0)  = 0 (B-11) 

The  tangential  velocity  of  a satellite  in  circular  orbit  is  given  by: 

V(,(0)  = (B-12) 


The  algorithm  described  in  Chapter  II  seeks  a value  for  the  initial 
mass  of  the  vehicle. 
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The  desired  final  conditions  are  determined  from  algebraic 


I 

/ 


relationships  as  in  Appendix  A.  The  final  eccentricity,  e£,  semi-latus 
rectum,  pf,  and  desired  transfer  angle  in  the  plane  of  the  final  orbit, 
0f,  are  given.  The  final  radius  is  therefore  given  as: 


^•f 


Pf 

1 + e£  cos 


(B-13) 


In  order  to  determine  Of  and  <j>f  in  the  initial  orbit  plane  we  must  use 
spherical  trigonometry.  The  spherical  triangle  which  describes  the 
relationship  between  the  initial  and  final  orbit  is  shown  in  Fig.  B-4. 


Fig.  B-4.  Spherical  Angles  for  Noncoplanar  Problem 
Napier's  rules  for  right  spherical  triangles  gives  the  following 
relationships: 


sin  = cos  (90-i)  cos  (180-^') 

sin  (#-90)  = tan  tan  (90-i) 


(B-14) 


\ 


Using  reduction  formulae  and  knowledge  of  </>q  we  can  obtain  the 
required  relationships: 


*#■1  = 90  - sin”^  (-sin  i cos  d'} 
H = cos“l  [tan  (90-<f>^)  cot  i] 


It  is  now  required  to  derive  expressions  for  the  velocities  at  the  final 
time.  The  expression  for  radial  velocity  shown  in  Appendix  A applies 
in  the  noncoplanar  case  as  well: 

, efhfSin»'f 


Where  hj  is  the  angular  momentum  of  the  final  orbit.  Relationships 
for  and  can  be  found  by  examining  expressions  for  angular 
momentum.  By  definition 


h = r X v 


which  in  spherical  coordinates  is: 


h = r V^e^  - r Vg  sin<#»e^ 


Another  expression  for  angular  momentum  in  the  terminal  orbit  is: 


h = h£  sin  i "E’x  + hf  cos  i e^ 


Eq  (B-19)  and  (B-20)  can  be  equated  after  the  coordinate  transforma- 
tion between  the  x-y-z  frame  and  r-^-0  frame  is  made.  The  transfor- 
mation matrix  given  by  three  successive  rotations  of  the  spherical 
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I 

i 

! 


1 


coordinate  system  through  angles  of  (90°-</>),  0,  and  90°  will  give  the 
appropriate  expressions  for  and  "e^  in  Eq  (B-ZO). 

= cos  9 sin  <t>  "ej.  + cos  9 cos<(>'e^  - sinO  eg  (B-21) 

"e^  = cos^'Sj.  - sin</>'e<^ 

The  expressions  that  result  from  substituting  Eq  (B-21)  into  Eq  (B-20) 
and  equating  that  to  Eq  (B-19)  is; 

Vg  = — (cos  i - sin  i cos  0£  cot  <(>^)  (B-22) 

^r 

sin  i sin  9 

Eqs  (B-13),  (B-15),  (B-l6),  (B-17),  and  (B-22)  are  the  expressions 
required  to  give  the  boundary  conditions  at  the  final  time. 
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Appendix  C 

Derivation  of  Necessary  Conditions 

The  presence  of  constraints  on  problem  variables,  i.  e.  mass 
and  time,  at  intermediate  points  (first  stage  cutoff,  staging,  and 
second  stage  ignition)  leads  to  a set  of  unique  necessary  conditions  for 
this  multipoint  boundary  value  problem  (MPBVP).  Those  conditions, 
as  derived  by  the  use  of  the  Calculus  of  Variations,  are  presented  in 
this  chapter  for  a complete  understanding  of  the  problem  formulation. 

Cost  Fimction 

The  problem  of  maximizing  the  payload  of  a space  vehicle  with 
a fixed  amount  of  fuel  and  structure  drop-off  weight  is  interpreted  as 
minimizing  the  negative  of  the  initial  mass.  With  that  interpretation 
the  performance  index,  or  cost  function,  as  described  in  Ref  5 is: 

J = - m(tQ)  (C-1) 

where  m is  the  mass  state  variable  and  is  the  initial  time.  The 
intermediate  boundary  conditions  that  must  be  met  in  the  problem 
appear  as  inequality  constraints  on  mass  and  thrusting  times  for  each 
of  the  two  stages.  Those  conditions  are: 

C - [m(ts")  - m(tg+)]^0  (C-2) 
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t,  - ti  > 0 

1 max  ^ — 


t2m.-oc  - <‘f  - '2>  i 0 


(C-3) 


where  C is  the  maiximum  allowable  change  in  mass  at  staging 

m(tg")  is  the  mass  just  prior  to  staging 

m(t  +)  is  the  mass  just  after  staging 
s 

ti  is  the  maximum  allowable  time  for  first  stage  burning 

1 rricix 

tj  is  the  first  stage  cutoff  time 

t is  the  maximum  allowable  time  for  second  stage  burning 

2 max 

is  the  second  stage  ignition  tinae 
tf  is  the  time  the  orbit  transfer  is  complete. 


The  end  conditions  which  must  be  met  are  the  desired  final  values  of 
the  states  (except  mass)  which  are  specified  a priori.  They  are 
denoted  by  the  vector; 


^ [x(tf)] 


(C-5) 


where  x represents  the  state  variables.  The  cost  function  is  minimized 
subject  to  the  differential  equations  in  the  state  variables: 


f(x,u,  t)  = X 


(C-6) 


where  u is  the  control  variable. 

An  augmented  cost  function  is  formed  by  adjoining  the  above 
mentioned  constraints  to  the  original  cost  function  by  the  use  of 
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lagrange  multipliers.  Due  to  the  possibility  of  discontinuities  in  the 
problem  the  differential  constraints  are  separated  into  four  intervals 
in  the  augmented  cost  function.  The  cost  function  is  now  stated  as: 


I 

j 


J = - m(to)  + i [x(t£)]  + )Tj  [C-m(tg  ) + m(tg''')  ] 


[•^Imax"*^!]  + ^^3  t‘^2max"*^f  + ^2^ 

r 

(fj^  - X)  dt  + A.T  (f^-x)  dt 

4 


(C-7) 


where  t9, 
the  state 


X,  and  TT  are  lagrange  multipliers,  and  f^,  i = I,  . , . 4 are 
equations  during  the  consecutive  phases  of  the  problem. 


The  optimization  problem  defined  by  this  performance  index  can 
be  reduced  to  an  MPBVP  by  applying  the  necessary  conditions  of  opti- 
mality which  will  be  derived  from  the  classical  Calculus  of  Variations. 


Necessary  Conditions 

The  necessary  conditions  are  derived  by  performing  a multi- 
variable  Taylor  series  expansion  of  J about  a nominal  trajectory. 
Higher  order  terms  are  neglected  and  the  first  variation  is  set  to  zero 
following  the  general  method  in  Ref  1 . Variations  of  the  form: 
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x(t)  = X*(t)  +8x(t) 


u(t)  = u*(t)  + 8u(t) 


(C-8) 


k(t)  = x*(t)  + 8x(t) 

■where  ( )*  indicates  a nominal  point,  give  the  following  result  after  a 
Taylor  Series  expansion  and  integration  by  parts  of  the  appropriate  terms: 

= [^x  (*^0)  dx(tf)  + 

-’^3!  dtf  + [ -^i-^m(t;)]  dm(r)  + 

+ [’"i  + dm(ts)  + [-7^2  + ‘^'^1 


+ ■ >’’^(tpdx(tp 

+ xT(t+)dx(t+)  + [xT(t;)f2  - xT(t+)f3]  dtg 
-x'^lt^)  dx  (ti)  + xT(t+)  dx  (t+) 


+J  ^ [(X  Tfj^  + A*^)  8x  + 

o 

+ f [ ®x  + ^2u  8u]  dt 

+ f ^ [(x'^f,  + x"^)  8x  + % «»i]dt 

•'i.T  -'X  ^ 


+ J^+  [(X^f4x  + »x  + 


(C-9) 
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Requiring  each  coefficient  of  the  nonzero  differentials  in  the  states 
(dx)  and  times  (dt)  to  be  zero  yield  the  necessary  conditions. 

Three  classes  of  trajectories  are  examined  in  this  thesis.  They 

are: 

Case  One--only  first  stage  fuel  is  completely  used. 

Case  Two--only  second  stage  fuel  is  completely  used. 

Case  Three- -both  first  and  second  stages  are  completely  used. 
The  necessary  conditions  will  first  be  developed  for  each  of  the 
specific  cases,  followed  by  a listing  of  the  conditions  common  to  all 
three  classes  of  trajectories. 

Case  One.  The  first  case  dictates  that  we  are  on  the  boundary  of 
the  first  stage  fuel  constraint  and  inside  the  boundary  region  of  the 
second  stage  constraint.  Examination  of  the  theory  of  constrained 
minimization  (Ref  1 ) provides  insight  regarding  the  lagrange  multi- 
pliers in  Eq  (C-7).  Being  on  the  first  stage  burn  constraint  indicates 
that  7>'2  nonzero  and  being  inside  the  second  stage  constraint  gives 
the  value  of  zero.  Burning  all  of  the  first  stage  fuel  also  indicates 
that  we  are  on  the  boundary  of  the  mass  constraint,  giving  a non- 
zero value.  We  also  know  that  a complete  burning  of  the  first  stage 
gives 

dm(t“)  = dm(t'*')  (C-10) 

s s 

since  the  structure  drop-off  weight  is  fixed.  The  resulting  necessary 
conditions  from  Eq  (C-9)  peculiar  to  case  one  are  therefore: 
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■JMlliLij'WiiiJIlwiMii.  ^ 


= 0 

v<‘;>  = ' ’i 

Hjitp  = H^(4) 


(C-11) 


where  for  i = 1,  4 

Case  Two.  For  this  class  of  trajectories  we  are  inside  the 
boundary  region  of  the  first  stage  burn  constraint  and  have  met  the 
constraint  for  the  second  stage  burning  period.  Constrained  minimiza 
tion  theory  dictates  that  and  are  zero  and  that  is  nonzero. 
Due  to  the  fact  that  the  first  stage  fuel  may  not  be  completely  used. 


dm(t”)  ^ dm(tg) 
s 

The  necessary  conditions  vmique  to  case  two  are: 


(C-12) 


H(tf) 

(f)  =A  (t+)  = 0 (C-13) 

m ' s'  m.'  s 

Hj(tj)  = H^itj) 

B^{t2)  - H3(t2)  = 

Case  Three.  The  third  case  indicates  that  we  are  on  the  boundary 
region  for  all  three  intermediate  constraints.  The  associated 
lagrange  multipliers  are  therefore  all  zero.  The  necessary 
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conditions  become; 


A (f)  =A  (r)  = 7T, 

m s'  m s'  1 


Hj(tj)  - H2(t^)  = -2 


H3(t2)  - H4(t2)  = - 


(C-14) 


General  Necessary  Conditions.  The  necessary  conditions  that 
must  be  met  for  this  problem  in  each  of  the  three  cases  can  be  deter- 
mined from  the  remaining  terms  in  Eq  (C-9).  Those  equations  are: 


A^(t)  continuous  for  all  x 


‘ <‘o)  " W'o>  ' + ‘ 


= xT(tf) 

X 


V(tf)  = 0 

d X 


(C-15) 


(C-16) 


(C-17) 


(C-18) 


(C-19) 


(C-20) 


(C-21) 


The  equations  that  correspond  to  Eq  (C-18)  are  commonly  referred  to  as 
the  transversality  conditions.  Those  equations  corresponding  to  Eq(C- 
20)yield  the  costate  equations  and  Eq  (C-21)  gives  the  optimality  condition 


on  the  control. 
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